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ARSDIGITA UNIVERSITY MONTH 0: 
MATHEMATICS FOR COMPUTER SCIENCE 



This course in concerned with the study of functions. For the first two weeks, we will 
study calculus, which will involve the study of some fairly complex functions. We will only 
skim the surface of much of this material. During the second two weeks, we will study 
linear algebra, which will involve the study of fairly simple functions - linear functions - in 
much greater detail. 

This is the outline of the course notes from September 2000. In the future, the syllabus 
may be slightly different, and so these notes may be out of date in a year's time. 

Lecture 1: Graphing functions. Exponential and logarithmic functions. 

1.1. What is a function? A function is a rule or a map from one set A to another set 
B. A function takes an element of the first set, the domain, and assigns to it an element of 
the second set, the range. We will use the following notation: 



In calculus, we will be interested in real- valued functions / : M ^ M. We can graph these 
functions. 

1.2. Graphing Functions. We can make a graph of a function y — f{x) by plotting the 
points {x, f{x)) in the xy-plane. 



We can make the graph roughly simply by choosing some values, plotting those vahics, 
and connecting the dots. We do have to be careful when we choose which x- values to plot. 



f,g,h,F,G are all functions. 



f : A ^ B is a function which maps x & A to f{x) G B. 




Figure 1. Here is how we graph a function y = 
f{x). We graph the points {x,f{x)). 
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1.3. Graphing tricks. If we know how to graph a function y = f{x), then we can build 
up a repertoire of functions by figuring out how this graph changes when we add constants 
to our equation: f{x) f{x + c), f{x) f{c ■ x), f{x) c ■ f{x), or f{x) f{x) + c. 
The graph changes by stretching and shifting, depending on the location of the constant. 
Tliis allows us to have a whole family of functions whose graphs we can draw, by having a 
few simple examples of functions in mind. 

1.4. Examples of Functions. The first examples are lines. To graph a line, we merely 
have to plot two points, and connect them with a straight line. The first example we have 
is the line y = x. A general line is of the form y = mx + b. The first constant m is called 
the slope and the second constant b is called the intercept. 

The next examples are polynomials. These are functions of the form y = anX^ + 
Qn-ix'^'^ + ■ ■ ■ + aix + Oo, where a„ 7^ 0. We say that the polynomial y = a^x"^ + a„_ia;"~^ + 
• • • + aix + ao has degree n. A line is a polynomial of degree less than or equal to 1. 

1.5. Exponentials and logarithms. Let a be some real number greater than 0. Then 
we can look at the exponential function y — and the logarithmic function y — log^ x. 

1.6. Trigonometric functions. There is a family of trigonometric functions, y = sin(a;), 
y — cos(a;), y = tan(a;), y = csc(a;), y = sec(a;), and y = cot(x). These are defined for 
angles x, measured in radians. These functions and their properties will be discussed in 
recitation. 

Lecture 2: Limits and continuity. Slopes and derivatives. 

This is the beginning of calculus in earnest. Single variable calculus addresses two major 
problems. 

1. Differential calculus addresses the question, what is the slope of the tangent fine to 

y — f{x) at the point P = {xq, /(xq))? 

2. Integral calculus, on the other hand, addresses the question, what is the area under 

the graph y = f{x) between x = a and x = bl 

This week, we will look at the first of these questions, and we will apply the answer to 
many real-world situations. 

2.1. Limits. Suppose we have a function y = f{x) defined on an interval around x = a, 
but not necessarily at a; = a itself. Suppose there is a real number L G M such that as x 
approaches a (from both the left and the right), f{x) approaches L. In this case, we say, 

lim f{x) = L or f{x) L as x ^ a. 
If no such number exists, then we say there is no limit. We will compute several examples. 

2.2. Continuity. We say that a function y — f{x) is continuous at x = a if 

1. /(a) is defined; 

2. lim^^^a exists; and 

3. hm^^a = f{o)- 

We will look at a couple of examples. 
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2.3. Slope and derivatives. The goal is to compute the slope of the tangent line to a 
point {x, f{x)). We do this by computing the slope of secant lines, and taking the limit as 
the secant lines approach the tangent lines. Thus, the derivative is 

dx Aa;^0 Ax Ax-'O Ax 

We will look at several examples of computing derivatives by this definition. 

2.4. Velocity and Rates of Change. In recitation, we will see how derivatives appear 
in the real world as the velocity of a position function, and as a rate of change. We will 
look at several problems of this nature. 

Lecture 3: Differentiation methods. 

In this lecture, we will see several methods for computing derivatives, so that we do not 
always have to use the definition of derivative in order to do computations. 

3.1. Rules for differentiation. Here is a list of differentiation rules. We will prove some 
of these in class, and use them to compute some difficult derivatives. 

1. = 0 for any constant c e R. 

2. ^(x") = n-x"-^ 

3. £{c ■ f{x)) = c ■ ^{f{x)), for any constant c. 
4- £{f{x)+9{x)) = i{f{x)) + i{g{x)). 

5. The product rule. 

^{f{x) ■ g{x)) = g{x) ■ ^{f{x)) + f{x) ■ ^{g{x)). 



6. The quotient rule. 



d ffix)\ g{x)-£{f{x))-f{x)-£{g{x)) 



7. The chain rule. 



where u = g{x). 
8. The power rule. 



dx \g{x)J (gix))^ 



d I n\ n-1 du 

dx dx 



9. ^(sin(x)) = cos(x). 

10. ^(cos(x)) = — sin(x). 

11. l(e^) = e^. 

12. £(a^) =ln(a) -a^. 

13. |(ln(x)) = i. 

14. £(log,(x)) = 

We will also discuss implicit functions and implicit differentiation. 
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3.2. Graphing functions. Derivatives can be used to graph functions. The first derivative 
teUs us whether a function is increasing or decreasing. The second derivative teUs us about 
the concavity of the function. This will be discussed in recitation. 

Lecture 4: Max-min problems. Taylor series. 
Today, we will be using derivatives to solve some real-world problems. 

4.1. Mctx-min problems. The idea with max-min problems is that a function f{x) defined 
on an interval [a, b] will take a maximum or minimum value in one of three places: 

1. where f'{x) — 0, 

2. at a, or 

3. at b. 

Thus, if we are given a problem such as the following, we know where to look for maximums, 
minima, the biggest, the smallest, the most, the least, etc. 

Example. Suppose we have 100 yards of fencing, and we want to enclose a rectangular 
garden with maximum area. What is the largest area we can enclose? <0 

We will discuss methods for solving problems like the following. 

4.2. Approximations. We can also use the derivative to compute values of functions. 
Suppose we want to compute \^A.02. Then we can take a linear approximation, the tangent 
line, of ^/x at a; = 4, and compute the approximate value of \/4.02 using the tangent line. 
Linear and quadratic approximations will be discussed. 

4.3. Newton's method. Often times, we would hke to know when a function takes the 
value 0. Sometimes, the answer to this question is obvious, and other times it is more 
difficult. If we know that f{x) has a zero between x = 1 and x = 2, then we can use 
Newton's method to determine the exact value of x for which f{x) = 0. This will be 
discussed in recitation. 

Lecture 5: Taylor Series. Antiderivatives. 
In this lecture, we will tie up loose ends of differential calculus, and begin integral calculus. 

5.1. Taylor series. The simplest functions that we have seen so far are polynomials. These 
are easy to evaluate and easy to differentiate. Suppose we have some more difficult function, 
y — f{x). Then we hope to approximate y — f{x) using a polynomial. Our first guess is 
the tangent line. We can add on higher terms using higher order derivatives. If we infinitely 
many terms, we end up with the Taylor series. The computation of the Taylor series will 
be discussed. 
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5.2. Antiderivatives. If we are given a function f{x), we say that another function F{x). 
is an antiderivative of f{x) if F'{x) — f{x). We will write 

J fix) dx^F{x). 

Notice that a constant can always be added to F{x)^ since = 0! We will compute several 
anti-derivaties. 

Taking antiderivatives can be quite difficult. There is no genearal product rule or chain 
rule, for example. We have some new rules, however. 



5.3. Change of variables. The change of variables formula is the antidifferentiation ver- 
sion of the chain rule. It can be stated as follows. 



j f{g{x))g'{x) dx^ j f{u) du. 



There is also an antidifferentiation version of the product rule. We will discuss this at a 
later time. 



Lecture 6: Area under a curve. Fundamental Theorem of Calculus. 

I claimed at the beginning of the course that integral calculus attempted to answer a ques- 
tion involoving the area under a curve. We will see how this is related to antidifferentiation 
today. 

6.1. Riemann sums. We can approximate the area between a function and the x-axis 
using Riemann sums. The idea is that we use rectangles to estimate the area under the 
curve. We sum up the area of the rectangles, and take the limit as the rectangles get 
more and more narrow. This leads us to the left-hand and right-hand Riemann sums. The 
left-hand Riemann sum is 



J a 



n-1 

f{x) dx = lim y f{xk) ■ Ax. 

Ao;— >0 ' 
fe=0 

The right-hand Riemann sum is 

/ f{x) dx = lim y^^fixk) ■ Ax. 

"'^ fe=l 

These both converge to the area we are looking for, so long as f{x) is a continuous function. 
We will make several computations using Riemann sums, and we will look at a few more 

examples of Riemann sums. 

These Riemann sums seem to have to relation to the antiderivatives discussed above, 
although we have used the same notation. It turns out that these two ideas are realted. 
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6.2. The Fundamental Theorem of Calculus. We will prove the following theorem. 

Theorem 6.2.1. If J{x) is a continuous function, and F{x) is an antiderivative of f{x), 
then 

I f{x) dx^F{h)-F{a). 

J a 

We will discuss the implications of this theorem, and look at several examples. Furthermore, 
we will discuss geometric versus algebraic area. 



Lecture 7: Second Fundamental Theorem of Calculus. Integration by 
PARTS. The area between two curves. 

7.1. Discontinuities and Integration. First, let's talk about one horrific example where 
everything seems to go wrong. This will be an example of why we need a function f[x) to 
be continuous in order to integrate. Consider the function 

V _ J 0 a; is a rational number, 
^ ' 1 1 x is an irrational number. 

We can't graph this function. The problem is, between every two rational points, there is 
an irrational point, and between every two irrational points, there is a rational point. It 
turns out that this function is also impossible to integrate. We will discuss why. 



7.2. The Second Fundamental Theorem of Calculus. When we were taking deriva- 
tives, we never really came across functions we couldn't differentiate. From our little bit 
of experience with taking indefinite integrals (the ones of the form J f{x) dx), it doesn't 
seem obvious that we can always integrate a function. The Second Fundamental Theorem 
of Calculus tells us that an integral, or antiderivative, always does exist, and we can write 
one down. 

Theorem 7.2.1. (Second FTC) If f{x) is a continuous function, and a is some real num- 
ber, then 

This is a by-product of the proof we gave yesterday of the First Fundamental Theorem 
of Calculus. 



7.3. Integration by Peirts. I talked about forming the "inverse" of the Product Rule 
from differentiation to make a rule about taking integrals. This is known as integration by 
parts, and the rule is 

Judv — u-v — Jv du. 

We will discuss this rule further, how to apply it, and some tricks for knowing which choices 
to make for u and v. 
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Figure 2. We're trying to compute the shaded area 
between f{x) and g{x) from x = a to x = b. 

7.4. The area between two curves. Finding the area between two curves looks pretty 
difficult. Suppose f{x) > g{x). See the figure. If we think about it, though, this is not a 
hard problem at all! When we think about it, all we have to do is take the area under f(x) 
on [a, b], and subtract from that the area under g{x) on [a, b]. Thus, the area between f{x) 
and g{x) on [a, b] is 

pb pb pb 

/ f{x) dx- g{x) dx= {f{x) - g{x)) dx. 

J a J a J a 

We will discuss some subtleties with this rule. 

Lecture 8: Volumes of solids of rotation. Integration techniques. 

8.1. Volumes of solids of rotation. Today, we'll talk about two methods to compute 
the volumes of rotational solids, volumes by discs and volumes by cylindrical shells. 

Let us consider for a moment the intuitive idea behind a definite integral. The definite 













c 


I X b 





integral f{x)dx is supposed to be an infinitesimal sum (denoted by J^) along the x- 
axis (denoted by dx) of slices with height f{x). The methods we will present to compute 
volumes of rotation will have a similar intuitive idea. 

8.LL Volume by slices. Consider a function f{x) from a to b. We can take the region 
between this curve and the a;-axis and rotate it around the x-axis. We get a solid of 
revolution. 
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The goal is to compute the volume of this solid. We have seen some examples of this 
already in recitation. The first method we will use to compute a volume like this is area by 




Figure 3. This shows the sHce at the point x. We 
say that this shce has area A{x). 

shces. The volume is the same as infinitesimally adding up the area A{x) of the slices, as 
shown in the above figure. So 

rb 

Volume = / A{x)dx . 

J a 

But this slice is just a circle with radius f{x), so A{x) — tt • = tt • f{xY. 

Volume = / nf{xYdx. 

J a 

8.1.2. Volume by shells. Again, we consider a solid of revolution. In this case, instead of 
slicing up our solid, we think of it as filled by cylinders, paper-thin cylinders. Thus, we will 
often use this method when the solid whose volume we wish to determine is described by a 
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function rotated around the y-axis rather than the x-axis. We can think of the volume as 
the infinitesimal sum of the surface area of the cylinders: 



Surface Area = 27irh 




Thus, when computing the volume by cylindrical shells, we have 



Volume 



2nxf{x) dx . 



Care is required to figure out what the appropriate height function is, and what bounds 
of integration to use. Furthermore, it is often useful to find a voume such as one of these 
by integrating with respect to y rather than x. We will see several examples of the cases 
to which these shces and cyhndrical shells apply. 

8.2. Integration using partial fractions. A rational function is the quotient of two 
polynomials, It is proper if 

deg(P(x)) < deg(Q(x)) 

and improper if 

deg{P{x)) > deg{Q{x)). 

An improper rational function can always be reduced to a polynomial plus a proper rational 
function. For example, 

(x^-l) + l 1 



x 



x^ 



x^ 



x^ 



, and 



{x' 



X 



x^ 



X' 



1 



X 



x^ 



Furthermore, we can "simplify" or combine terms: 



+ 



6 



2(a; - 2) + 6(a; + 3) 8x + U 



X + 3 X — 2 



(a;-2)(a; + 3) 



x^ + X — 6 



How can we go backwards? This is most easily illustrated with an example. If we start 
with 

2x-n _ A B 
{x + 2){x - 3) ~ x + 2 ^ x-3 ' 
how do we find A and B7 We know that A and B must satisfy 

A{x -3) + B{x -2) = 2,7; -11, so 
{A + B)x + {2B - 3A) = 2a; - 11 , 
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so 



A + B = 2, 
-3A + 2B = -11. 

We can solve this system, or we can substitute x = 3,—2 in A{x — 3) + B{x + 2) to get 

^(3 -3) + 5(3 + 2) = 2-3-11 
5B = -5 
B = -1, 

and 

A{-2-3) + B{-2 + 2) = 2(-2)-ll 
-5A = -15 
A = 3. 

So 

2a; - 11 3 1 



{x + 2){x-3) x + 2 x-3' 
What is this good for? Suppose we want to integrate 

2x-n 



I 



dx . 



{x + 2){x-3) 

We might try to substitute u = x^ + x — 6, but we would soon see that {2x — 11) dx is not 
du. So instead, we can split this up: 



2a;- 11 



{x + 2){x - 3) 



dx 



dx 



x + 2 
u = x + 2 



dx 



X — 3 
w = X — 3 



du 



dx dw — dx 

dw 

31n b + 21 — In la; — 31 + c 



= [ — du — [ — I 
J u J w 




A couple of cautionary notes about partial fractions: 

1. If there is a term (a; — a)"* in the denominator, we will get a sum of the form 

Am, 



Ai _^ A2 



+ ■■■ + 



a; — a 



X — a \x — a)^ 

where the A-i are all constants. 
2. If there is a quadratic factor ax^ + 6a; + c in the denominator, it gives a term 

Ax + b 
ax"^ + bx + c 



©Tara S. Holm, MMI. May be distributed as per http://opencontent.org. 



Outline 



11 



8.3. Differential Equations. Finally, if there is time, we will talk a bit about difTerential 
equations. A differential equation is an equation that relates a function and its derivatives. 
For example, consider 

dy 

— = xy. 

dx 

The goal is to determine what ?/ is as a function of x. If we treat the term ^ as a fraction, 
we can separate the variables and put all terms with y and dy on the left-hand side, and 
all terms with x and dx on the right-hand side. Upon separating, we get 



— — X dx. 



y 

We can now integrate each side to get 



dy f x^ 
In ||/| = / — = X dx = — + c. 



y J 2 

Thus, by raising each side to the power of e, we get 

y = C-e^\ 

where C is a constant. Differential equations that can be solved in this manner are called 
differential equations with variables separable. There are many other differential equations 
which can not be solved in this way. Many courses in mathematics are devoted to solving 
different kinds of differential equations. We will not spend any more time on them in this 
course. 



Lecture 9: Trigonometric substitution. Multivariable functions. 

9.1. Trigonometric substitutions. Today, we'll discuss trigonometric substitutions a 
little further. To use the method of trigonometric substitution, we will make use of the 
following identities. 

l-sin^^ = cos^^, 
1 + cot^^ = csc^^, 
sec^ e-1 = tan^ 9 . 

If we see an integral with the function \/ — x^ in the integrand, with a a constant, we 
can make the substitution x — a%\TiQ . Then we have 

a? — x'^ — \/ — sin^ Q — V a? cos^ 9 — a cos 9. 

Notice that in the case, dx = a cos 9 d9. We make similar substitutions in the remaining 
cases. If we see a function of the form \J x"^ — in the integrand, we can make the substi- 
tution X = a sec 9. If wc sec a function of the form \/a?M-c? in the integrand, we can make 
the substitution x = a cot 9. 
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9.2. Cavalieri's Principle. Next, we'll discuss Cavalieri's Principle. This relates to the 
disc method we discussed yesterday, although Cavalieri, a contemporary of Gahleo, stated 
the result centuries before the advent of calculus. We will see that this principle is true 
using calculus. 

Theorem 9.2.1. If two solids A and B have the same height, and if cross- sections parallel 
to their bases have the same corresponding areas, then A and B have the same volume. 

This is clearly true, if we take the volumes by discs. Since the areas of the discs are the 
same, the volume that we compte when we integrate will be the same. 

9.3. Multivariable Functions. Finally, there will be a brief introduction to multivariable 
functions. We will discuss graphing multivariable functions and multiple integrals, as time 
permits. 

Given a multivariable function z — f{x,y), we can graphs its level curves. These are 
the curves that we obtain by setting z = a, for constants a. As a varies, we get all of the 
level curves of the function, which we graph in the xy-pla.ne. For example, we can draw 
the level curves of the function z — +y^. These will just be circles a — x^ +y^. Using 



the level curves, we can get a three dimensional view of the function. Now instead of a 



tangent line to a curve, we have a tangent plane to a surface. We can still compute partial 
derivatives. But there are many direections in which we can differentiate. We now have 
directional derivatives. Finally, we can also compute the volume under a surface z = f{x, y) 
over a region R. To do this, we can break R up into many small squares, and add up lots 
of rectangular solids, taking their height to be some avearge value of the function. We can 




Figure 4. This shows the level curves of the func- 
tion z = + y'^. 




Figure 5. This shows the three dimensional graph 
of the function z = + y^. This graph is known as 
a paraboloid. 



©Tara S. Holm, MMI. May be distributed as per http://opencontent.org. 



Outline 



13 



make exactly the same Riemann sums as before, but this time instead of just dx, we have 
dA — dx dy — dy dx. 

Lecture 10: Linear algebra 

Today, we make a total break from calculus and start linear algebra. We (probably) 
won't take another derivative or integral this month! In linear algebra, we will discuss a 
special class of functions, linear functions 

/ : ^ M™. 

We will begin by discussing vectors, properties of vectors, and matrices. 

10.1. Vectors. We have vectors in W^. For example, 

1 
2 
1 
0 
3 

In general, we will think of vectors in as column vectors, but we may write them as row 
vectors in order to save space. We will use v and w to represent vectors. We add vectors 
component-wise. That is. 



" 1 " 




2 " 




" 1 + 2 " 




3 " 


2 




0 




2 + 0 




2 


1 


+ 


5 




1 + 5 




6 


0 




-4 




0-4 




-4 


3 




1 




3 + 1 




4 



We can multiply a vector by a real number a & R, a scalar, and we do this component- wise 
as well. We cannot multiply two vectors, but we can take their dot product. 







' yi ' 






1/2 






. Un . 



E 

1=1 



Xi ' yi 



The dot product of two vectors is a real number. The length of a vector f is | |f | 
A unit vector is a vector v with ||f || = 1. If we have two vectors v and w in W' 
angle 9 between them can be computed by 

V ■ tr 

COSd = 



= \/v ■ V. 
then the 



\\v\\ ■ \\w\\ 

This is known as the Law of Cosines. Notice that if the vectors are perpindicular, then 
the cosine is 0, so the dot product of the vectors must be 0. We also have the Schwartz 
Inequality, which says 

V ■ w < \ \v\ \ ■ \ \w\\. 
Note that this follows from the Law of Cosines. 
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10.2. Matrices. A matrix is an m x n array. An mxn matrix has m rows and n columns. 
We will use a matrix to store data. The entry of a matrix is the entry in the i*^ rown 
and j*^ column. Consider the matrix 

" 342 216 47 " 
312 287 29 ' 

The first row might represent sales in August 2000, and the second sales in August 1999. 
The first column might reperesent gallons of milk, the second loaves of bread, and the third 
heads of lettuce. So long as we remember what each row and column represents, we only 
need keep track of the entries in each position. 

As there were vector operations, we also have matrix operations. We can add two ma- 
trices, so long as they are the same size. As before, addition is entry-wise. For example, 

1-2 2 + 3 " 
-2 + 7 0 + 4 ■ 

We can also multiply a matrix by a scalar. Again, this is entry-wise. Finally, we can 
multiply two matrices together, so long as they have the right sizes. We can multiply an 
mxn matrix by an n x p matrix to get a m x p matrix. The way this happens is if we are 
multiplying A - B to get C, then the {i,j) entry of C is the dot product of the i^^ row of A 
with the j*^ column of B. Notice that the dimension restrictions on A and B are such that 
these two vectors are the same size, so their dot product is well-defined. For example. 



12" 




"-23" 




' l-(-2) + 2-7 -2 -(-2) + 0-7 


-2 0 




7 4 




1-3 + 2-4 -2 •3 + 0-4 



We can also use a matrix to represent the coefficients of linear equations. For example, 
if we have the system of equations 

f 342a; + 216y + 47z = 1041 
\ 312x + 2877/ + 29^ = 998 ' 

we can find values of x, y and z such that both equations are true. If we let a; = 2, y = 1 
and z = 3, the equations are true. We will use matrices to find all possible solutions of 
these equations. We can let a matrix hold the coefficients of these equations. We can write 
down this system of equations using the matrix eqation 

" 342 216 47 
312 287 29 

We will often refer to this form as A ■ x = b, where A is the coefficient matrix, x is the 
vector of variables, and b is the desired solution vector. 

Lecture 11: Gaussian elimination. Matrix operations. Inverses. 

Today, we will further explore matrices. 
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X 

y 

z 



1041 
998 
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11.1. Gaussian Elimination. First, we will learn an algorithm for solving a system of 
linear equations. This Gaussian elimination will come up later in the algorithms course, 
when you have to figure out its time complexity. Suppose you have the following system of 
equations. 

( X + 2y - z = 2 
\^ X + y + z — 3 ' 

Subtracting the first equation from the second gives a new system with the same solutions. 

( X + 2y - z ^ 2 
\ -y + 2z ^ 1 ■ 

This system has infinitely many solutions. We can let the variable z take any value, and 
then X and y are determined. The existence of infinitely many solutions can also occur in 
the case of three equations in three unknowns. If a coefficient matrix A is n x n, or square, 
then we say that it is singular ii A ■ x = b has infinitely many solutions or no solutions. 

Before we give the algorithm for Gaussian elimination, we will define the elementary row 
operations for a matrix. An elementary row operation on an m x n matrix A is one of the 
following. 

1. Swapping two rows. 

2. Multiplying a row by a constant c G M. 

3. Adding c times row i of A to ci times row j oi A. 

Now, given a system of equations such as 

( X + 2y - z ^ 2 
\ X + y + 2 = 3' 

the augmented matrix associated to this system is 



" 1 


2 


-1 


2 " 


1 


1 


1 


3 



In general, for a system A - x — b, the augmented matrix is [ A | 6 ] . We say that a matrix 
is in row echelon form when it satisfies the following conditions. 

1. All rows consisting entirely of zeroes, if there are any, are the bottom rows of the 
matrix. 

2. If rows i and i + 1 are successive rows that do not consist entirely of zeroes, the the 
leading non-zero entry in row i + 1 is to the right of the leading non-zero entry in row 
i. 

In Gaussian elimination, we perform elementary row operations on the augmented matrix. 
We perform the operations to both sides of the matrix. We write the vertical bar simply to 

remind ourselves that the matrix is augmented. The goal of Gaussian elimination is to put 
the augmented matrix into row echelon form. This then allows us to solve the associated 
system of equations. We will discuss the precise details in class. 
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11.2. Inverse matrices. Finally, we will compute the inverse of a matrix. Given a square 
matrix A, its inverse is a matrix A~^ such that 

A- A-^ = A-^ ■ A = /, 

where / is the identity matrix. To compute the inverse, we take the augmented matrix 
[ A I / ] , and perform elementary row operations on it until it is in the form [ / | i? ] . 
Then the matrix B is the inverse of A. We will talk about how to prove this in class. 

A square matrix is not always invertible. It turns out that a matrix is not invertible if 
and only if it is singular. 

Lecture 12: Matrices. Factorization A = L - U . Determinants. 

12.1. Facts about matrices. Here are some interesting facts about matrices, their mul- 
tiplication, addition, and inverses. Each of the following is true when the dimensions of the 
matrices are such that the appropriate multiplications and additions make sense. 

1. It is not always true that AB — BA\ 

2. There is an associative law, A{BC) = {AB)C. 

3. There is also a distributive law A{B + C) = AB + AC. 

4. (AB)-^ = B-^A-\ 

There are also some important types of matrices. We say that an entry of a matrix is on 
the diagonal if it is the (i, i) entry for some i. A matrix is a diagonal matrix if all of its 
non-diagonal entries are zero. For example, the matrix 



r d^ 


0 


0 . 


.. 0 


0 




0 . 


.. 0 


0 


0 




.. 0 


. 0 


0 


0 . 


■ ■ dn 



is a diagonal matrix. Diagonal matrices have some nice properties. They commute with 
each other (that is, CD = DC), and they are easy to multiply. The identity matrix / is a 
special diagonal matrix: all its diagonal entries are 1. Another type of matrix is an upper 
triangular matrix. This is a matrix which has zeroes in every entry below and to the left 
of the diagonal. A lower triangular matrix, is a matrix which has zeroes in every entry 
above and to the right of the diagonal. Finally, the transpose of an m x n matrix A is an 
n X m matrix A^, where the rows and columns are interchanged. In this operation, the 
(i,_7)-entry becomes the (j, i)-entry. If a matrix has A — A^, then we say A is symmetric. 
Notice that (^1^)^ = A, and that {ABf = S^^t 

12.2. Factorization of A. We will discuss the factorization 

A^L-U 

of a matrix A into a lower triangular matrix L times an upper triangular matris U. It turns 
out that we already have the matrix U. This is the matrix that is the result of Gaussian 
elimination. If we keep track of our row operations, then we can also determine L. Note 
that if you need to perform a row exchange during elimination, then the factorization will 
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be of the form P ■ A — L ■ U. We will discuss how to find L and P if necessary in class. 
Finally, in the case that ^4 is a symmetric matrix, the factorization A — L-U becomes 

A^L-D-L^, 

where D is a diagonal matrix, and L is lower triangular (and hence is upper triangular). 
Our ultimate goal is to factor a matrix A into the form 

A^C-D-C-\ 

where C is an invertible matrix, and D is a diagonal matrix. This factorization will allow 
us to compute powers of A easily, since 

12.3. Determinants. As the name suggests, the determinant of a matrix will determine 
something about the matrix: if det{A) — 0, then A is a singular matrix. If A is a 2 x 2 
matrix, then 

a b 
c d 

and det(A) = ad — be. We will define the determinant for a larger square matrix recursively. 
Some important properties of the determinant are the following. 

1. det(/) = 1 for / the identity matrix. 

2. If two rows of A are the same, the det(A) = 0. 

3. If a row of A is all zeroes, then det{A) = 0. 

4. If is a diagonal matrix, then det{D) is the product of the diagonal entries. 

5. If A is upper or lower triangular, then det{A) is the product of the diagonal entries. 

6. det{AB) = det{A)det{B). 

Lecture 13: Vector spaces. 

13.1. Application of factorization. We have seen three algorithms thus far. They can 
be described as follows. 

1. Gaussian elimination. 

(a) Ax — b gives an augmented matrix [ A | 6 ] . 

(b) Use row operations to transform [ >1 | 6 ] into [ | c ] , where U is an upper 
triangular matrix. 

(c) Back solve to find solution(s) x. 

2. Gauss- Jordan elimination. 

(a) Ax — b gives an augmented matrix [ A | 6 ] . 

(b) Use row operations to transform [ A | 6 ] into [ / | o? ] . 

(c) Read off answer (s) x = d. 

3. A — L-U factorization. 

(a) Ax — b. Just consider A for a minute. 

(b) Factor A = L-U. 

(c) Ax = L- U{x) = b. 

(i) Forward substitute Lz = b. 

(ii) Back solve Ux — z. 
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We will discuss each of these algorithms, and see why the last algorithm is the fastest. 

13.2. Abstract vector spaces. We are not going to take a step towards the more abstract 
side of linear algebra. We are going to define a vector space, something like M". A vector 
space (over R) is a set V along with two operations, + addition of vectors and • scalar 
multiplication, such that for all vectors v,w, z &V and scalars c, d e M, 

1. V + w ^ V and v -\- w = w + v. 

2. {v + w) + z = V + {w + z). 

3. There is a zero-vector 0 such that 0 + v = v. 

4. For every v, there is an additive inverse —v such that v + {—v) — 0. 

5. c ■ V E V and c • {d • v) — {c ■ d) • v. 

6. c ■ {v + w) = c ■ V + c ■ w. 

7. {c + d)-v = c- v + c- w. 

Of these, the most important ones are numbers 1 and 5. Most of the other properties 
are implied by these two. The most important thing to check is that for any two vectors 
v,w E V, any linear combination c - v + d ■ w of them is also in V. 

In class, we will see many examples of vector spaces. We will continue our discussion 
of vector spaces by defining subspaces. A subset W Q V of a. vector space V is a vector 
suhspace if it is also a vector space in its own right. A map L : Vi — > V2 is a linear 
transformation if L{c • v + d ■ w) — c ■ L{v) + d ■ L{w). The kernel or nullspace of a linear 
transformation is the subset Wi C Vi such that L{w) = 0 for all w G Wi. We write the 
kernel as ker(L) or N{L). The nullspace of a vector space is a subspace. 

Lecture 14: Solving a system of linear equations. 

14.1. More vector spaces. First, we will go over several examples of vector spaces. The 
most important are the row space and column space of a matrix. Given a matrix A, the 

row space R{A) is the vector space consisting of all linear combinations of the row vectors 
of A. The column space C{A) is the vector space consisting of all linear combinations of 
column vcctos of A. 

14.2. Solving systems of lineeir equations. We will discuss the complete solution to a 
systme of linear equations, 

A- X ^b. 

To find the complete solution, we must put the matrix into reduced row echelon form. A 
matrix is in reduced row echelon form if it satisfies the following conditions. 

1. All rows consisting entirely of zeroes, if there are any, are the bottom rows of the 
matrix. 

2. If rows i and i + 1 are successive rows that do not consist entirely of zeroes, the leading 
non-zero entry in row i -|- 1 is to the right of the leading non-zero entry in row i. 

3. If row i does not consist entirely of zeroes, then its leading non-zero entry is 1. 

4. If a column has a leading non-zero entry (a 1!) in it in some row, then every other 
entry in that column is zero. 

The columns that contain a leading entry are called pivot colurrms. The columns that do 
not contain leading entries are called free columns. The rows that arc not entirely zeroes 
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are pivot rows. The rank of a matrix A is the number of pivot columns (or rows) that it has 
when it is transformed into reduced row echelon form. We write rank(A) = r. The pivot 
variables are the variables that correspond to the pivot columns. The free variables are the 
variables corresponding to the free columns. 

To find a particluar solution Xp to Ax = b, we put [ A | 6 ] into reduced row echelon 
form. Then we set the free variables equal to zero, and read off the values of the pivot 
variables. This is one solution. If a system has free variables, then it has infinitely many 
solutions. Every solution is of the form Xp + x^, where a:„ is a vector in the nuUspace of 
A, that is Axn = 0. It turns out that we can write down a general solution using the free 
variables. We will discuss this precisely in class. 

14.3. The fundamental subspaces. We will discuss the four subspaces associatied to an 
m X n matrix A. The four subspaces are the row space R{A), the column space C{A), the 

rp 

nuUspace N{A), and the left nuUspace N{A^). 

Lecture 15: Linear independence. Bases. 

15.1. Independence and dependence. Spanning sets. Let y be a vector space. Sup- 
pose vi, ■ ■ ■ ,Vn & V is some finite set of distinct non-zero vectors. A linear combination of 
these is a sum 

ClVi + C2V2 H h CnVn, 

where the q are real numbers. We say that the set of vectors {fi,--- is linearly 
independent if no non-zero comibination of the vectors is zero. That is, if 

ClVi + C2V2 H h CnVn = 0 

for real numbers Cj, then all the Cj must be zero if the Vi are independent. 

We say that the set of vectors {vi, • • • , Vn} is linearly dependent if there is some non-zero 
linear combination of them that is zero. This means that one of the vectors is a linear 
combination of some of the other vectors in the subset. That vector depends on the other 
vectors. 

The rank of a matrix can be viewed using the idea of independence. The rank is the 
number of linearly independent rows (or columns) of a matrix. 

We say that a subset of vectors {vi, • ■ ■ ,Vn} in V spans V if every vector v E V is a, 
linear combination of these vectors. We will also call this subset a spanning set. 

15.2. Bases and dimension. A basis is a set of vectors {vi, • • • ,Vn} in V that is both 
linearly independent and is a spanning set. 

Theorem 15.2.1. // {vi, • • • , Vn} and {wi, • • • , Wm} are both bases for a vector space V , 
then n — m. 

We will prove this theorem in class. This allows us to define dimension. The dimension 
of a vector space V is the number of vetors in a basis. It turns out that the dimension 
of the vector space IR" is n, so this definition of dimension conincides with the notion we 
already have. We will prove this in class. 

We now return our attention to the fundamental subspaces of a matrix. Suppose that 
A is an m X n matrix with rank r. Then it turns out that the dimension of the four 
fundamental spaces are as follows. 
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1. The row space R{A) C M"^ has dimension r. 

2. The column space C{A) C has dimension r. 

3. The nullspace N{A) C has dimension n — r. 

4. The left nullspace N[A^) C M™ has dimension m — r. 

Lecture 16: Orthogonality and orthonormality. Gram-Schmidt 

algorithm. 

16.1. Orthogonality. We say that a pair of vectors v and w in R" is orthogonal or per- 
pindicular if their dot product is zero, v ■ w — Q. In vector spaces that are not R", we must 
define what we mean by dot product before we can define orthogonal. Thus, we will restrict 
our attention to V + M". 

A basis B = {vi, . . . , Vn} of R" is an orthogonal basis if every pair of vectors is orthogonal. 
That is, Vi ■ vj — 0 whenever i j. For example, the standard basis of R^ is the basis 




This is an orthogonal basis. 

Finally, two subspaces Wi and W2 inside V are orthogonal subspaces if for every wi e Wi 

and W2 G W2, Wi and W2 are orthogonal. Suppose that A is an m x n matrix with rank 
r. Then the row space R{A) and the nullspace N{A) are orthogonal subspaces of M™, and 
the column space C{A) and the left nullspace A^(y4^ arc orthogonal subspaces of R". 

16.2. Orthonormality. We say that a basis B = of R" is an orthonormal 

basis if every pair of vectors is orthogonal, and if each vector has length 1. That is. 

The standard basis for R^ described above is in fact an orthonormal basis. Orthonormal 
bases are particularly easy to work with because of the following theorem. 

Theorem 16.2.1. If B — {qi, . . . , g„} is an orthonormal basis of R" and v & MT' is any 

vector, then 

V = CiQi H h CnQn, 

where Ci = v ■ is the dot product of v and q^ for all i. 

We say that a matrix is an orthogonal matrix if its columns form an orthonormal basis 
of the column space of the matrix. We will be particularly interested in square orthogonal 
matrices. The name orthogonal is a bit of a misnomer. Unfortunately, the name has 
mathematical roots, and so although we would like to call this an orthonormal matrix, we 
will not. Suppose 





T 


T 


T 


Q = 




q2 ■ 


• (In 






L 





is an orthogonal square matrix. Then Q is invertible, and furthermore, Q ^ = ■ 
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16.3. Gram-Schmidt algorithm. Suppose B = {vi, . . . ,Vn} is any basis of M". We 
will use the Gram-Schmidt algorithm to turn this into an orthonormal basis {qi, . . . , g„} 

such that the subspace spanned by {vi, . . . ,Vk} is the same as the subspace spanned by 
{qi, . . . ,qk} for aU k. The algorithm proceeds as follows. First we will turn B into an 
orthogonal basis. We let 

Wi = Vi. 

Next, the second vector W2 should be perpindicular to wi. We define it to be 



W2 



V2 ■ Wi 

V2 Wi. 

Wi ■ Wi 



We continue in this manner, and let 

w. 



Vi ■ Wi 
Vi Wi 



Vi ■ W2 



-W2 



Vi ■ Wi_i 



Wi ■ Wi W2 ■ y'2 Wj-l ■ Wi^i 

At this stage, we could stop, and we would have an orthogonal basis {wi, . . . , Wn} such that 
the subspace spanned by {vi, . . . , Vk\ is the same as the subspace spanned by {wi, . . . , Wk\ 
for all k. To turn this into an orthonormal bases, we normalize, and set 

Wi 

This basis {qi, . . . , g„} is the desired orthonormal basis. We will prove in class that this is 
indeed an orthonormal basis. 

16.4. Change of basis matrices. Let B be the standard basis of M". That is, B = 
{ei, . . . , e„}, where Cj is the vector that has a 1 in the i*^ position, and zeroes elsewhere. 
Let B — {vi, . . . , Vn} and B — {wi, . . . , w„} be two other bases of R". We will now discuss 
the change of basis matrices. Suppose we have a vector 



V — 



Cl 



Then we can write this vector as 

Cl 



CiCi H h CnCn = / • 



Cl 



We may want to write our vector in terms of B. To do this, we need to find constants 
bi, . . . ,bn such that 











T 


T 


= bivi + ■ 


■ ■ + bnVn = 


Vi ■ 
















we g 


et that 










■ fci " 




' T 


T " 


-1 








Vi ■ ■ 












L 









■ &1 ■ 







Cl 
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Thus, we say that the change of basis from B to B is 





T 


T 




Vi ■ 






± 





The change of basis matrix from S to S is 





T 


T 










± 





Finally, if we want to change basis from B to B, then the change of basis matrix is 

Wc seem to have multiplied in the wrong order. We will discuss this in class, and compute 
several examples. 

Lecture 17: Eigenvalues, eigenvectors and diagonalization. 

17.1. Eigenvalues and eigenvectors. For a matrix A, the eigenvalues of A are the real 
numbers A such that 

A • V = \ • V 

for some vector v. Let's explore this equation further. 

A-v ^ X-v = {X- 1) -v 

Av-XIv = 0 
{A-XI)-v = 0 

This last line will have a nonzero solution v if and only if 

det(^ - XI) = 0. 

Thus, we find the eigenvalues of A by finding roots of the polynomial 

det(^ - XI) = 0. 
An eigenvector of a matrix A is a vector v such that 

A ■ V — X ■ V, 

for some real number A. Once we know what the eigenvalues are, we can find associated 
eigenvectors. 

Eigenvalues do not necessarily exist, other than 0. They are the roots of a degree n poly- 
nomial, and so there are at most n eigenvalues, and there can be no non-zero eigenvalues. 
If wc have n eigenvalues, they may not be distinct. If we have n distinct eigenvalues, then 
if we choose one eigenvector for each one, these will form a basis for MJ^. 
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17.2. Diagonalization. If A is a matrix, and there is a diagonal matrix D and an invertible 
matrix P such that 

A^P-D-P-^, 

then we say that A is diagonalizable. 

Theorem 17.2.1. A matrix A is diagonalizable if it has n distinct eigenvalues. In this 
case, if Xi, . . . , Xn are the eigenvalues and vi, - ■ ■ ,Vn are the eigenvectors, then 

Ai 0 0 - - 



T 



T 

Vn 
± 



0 
0 



■•• 0 

0 A„. 



T 



T 

Vn 
± 



In the case when the eigenvalues are not distinct, the matrix may or may not be diago- 
nalizable. In this case, there is still a canonical form known as the Jordan canonical form. 
We will discuss this in class. Some matrices may be diagonalizable over M, but they are 
over C. We will see an application of diagonalization tomorrow. 



Lecture 18: Fibonacci numbers. Markov matrices. 

18.1. Fibonacci numbers. In the year 1202, the Italian mathematician Fibonacci pub- 
hshed a book in which he defined a certain sequence of numbers. The so-called Fibonacci 
numbers are defined by the relation 

Fn — Fn-i + Fn-2, 

and Fq — 1 and Fi = 1. This sequence is 

0,1,1,2,3,5,8,13,21,34,55,.... 

We can use diagonalization to write a closed form equation for the n}^ Fibonacci number. 
If we put two values in a vector, 

Fn-1 
Fn-2 



then we can write 



Thus, 



1 1 
1 0 



Fn-1 




' F^_ 


1 + Fn-2 




Fn-2 _ 






Fn-1 






' 1 1 " 


n 


' F, ' 




Fn+1 


1 0 




L ^0 




Fn 





Fn-1 



So if we could easily compute 



1 1 
1 0 



we could easily compute the n Fibonacci number. If we diagonalize the matrix 

" 1 1 " 
1 0 ' 
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then we can easily compute its n* power. We will do this in class, and come up with the 
formula 



18.2. Markov matrices. We will briefly discuss Markov matrices. This is a place where 
linear algebra can be applied. A Markov matrix represents a scenario where certain outcome 
occur with various probabilities, and these probabilities change over time. In some cases, 
a steady-state is approached, and we will see exactly when. We will use eigenvectors and 
eigenvalues to determine what that steady-state is. 
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